QUANTUM INTEGRALS OF MOTION FOR VARIABLE QUADRATIC 

HAMILTONIANS 
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' Abstract. We construct integrals of motion for several models of the quantum damped oscillators 

in a framework of a general approach to the time-dependent Schrodinger equation with variable 
^ ' quadratic Hamiltonians. An extension of the Lewis-Riesenfeld dynamical invariant is given. The 

time-evolution of the expectation values of the energy related positive operators is determined for 
the oscillators under consideration. A proof of uniqueness of the corresponding Cauchy initial value 
, problem is discussed as an application. 



Ph. 



1. An Introduction 



Evolution of a nonrelativistic quantum system from a given initial state to the final state is 
governed by the (time-dependent) Schrodinger equation. Unfortunately, its explicit solutions are 
available only for the simplest Hamiltonians and, in general, one has to rely on a variety of approx- 
^ , imation, asymptotic and numerical methods. Luckily among the integrable cases are the so-called 
O quadratic Hamiltonians that attracted substantial attention over the years in view of their great 
^ ■ importance to many advanced quantum problems. Examples can be found in quantum and phys- 
t:J- • ical optics [21], [Zl], [TTi] , pTB] . physics of lasers and masers [Hg] . [112], [TM] , [Tig] , molecular 
', spectroscopy [H], quantum chemistry, quantization of mechanical systems [21], [IS], [IS], [IZ], [5U] . 
^ ; [75], [77] and Hamiltonian cosmology [H], [H], [S2], [SB], [HI], [IE], [HI], [I2S], [I2S]. They include 
g ; coherent states [HZ], [H] and Berry's phase [7], [S], [IB], [SZ], [Bl], [IDS], asymptotic and 

^ ■ numerical methods [51], [68], [78], |103] . |107] . charged particle traps [91] and motion in uniform 
• ^ '. magnetic fields [26], [29], [39], [80], [88], [89], [91], [97], polyatomic molecules in varying external 
^ ■ fields, crystals through which an electron is passing and exciting the oscillator modes and other 
interactions of the modes with external fields ^0]. Quadratic Hamiltonians have particular applica- 
tions in quantum electrodynamics because the electromagnetic field can be represented as a set of 
forced harmonic oscillators [TU], [SO], [3S], [S3], [SS], and |1U2] . Nonlinear oscillators play a central 
role in the novel theory of Bose-Einstein condensation [5D] based on the nonlinear Schrodinger (or 
Gross-Pitaevskii) equation [HS], [HZ], [Z3], |119] . 

The one-dimensional Schrodinger equation with variable quadratic Hamiltonians of the form 
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where a{t) , b (t) , c{t) , and d {t) are real- valued functions of time t only, can be integrated in the 
following manner (see, for example, |2S], EZ], [2B], [M], iH, [S3], [M], HM], [ISS], USE], HSZ], 
|149] . and |150] for a general approach and some elementary solutions). The Green functions, or 
Feynman's propagators, are given by [26], |136] : 

V; = G'(x,y,t) = ^. g<aWx2+/3Wx,+7W,^)^ ^1^2) 
y 27rz/i (t) 



where 



a{t)= . - TTTT^^ (1-3) 



4a (t) ^1 (t) 2a (t) ' 
hit) 



P (t) = h (t) = exp ( - / (c (r) - 2d (r)) rfr ) , (1.4) 







fiWA'W 2o(0) 7„ ^ 

and the function /i (t) satisfies the so-called characteristic equation 

fx" - T (t) fi' + Aa (t) fi = (1.6) 

with 

r{t) = — -2c + Ad, a{t) = ab- cd + d"^ + ^ (- (1.7) 
subject to the initial data 

/i(0)=0, (0) = 2a (0) ^ 0. (1.8) 

(More details can be found in Refs. [26], |136] and a Hamiltonian structure is considered in Refs. [7], 
[28].) Then, by the superposition principle, solution of the Cauchy initial value problem can be 
presented in an integral form 

/oo 
G{x,y,t) ip{y) dy, Mm {x,t) = ip (x) (1.9) 
-oo *^o+ 

for a suitable initial function on M (a rigorous proof is given in Ref. |136] and uniqueness is 
analyzed in this paper). 

We discuss integrals of motion for several particular models of the damped and generalized 
quantum oscillators. The simple harmonic oscillator is of interest in many quantum problems 
|5U] . [50] . |1U2] . and |129] . The forced harmonic oscillator was originally considered by Richard 
Feynman in his path integrals approach to the nonrelativistic quantum mechanics [16], [17], [H], 
[19], and [50]; see also [93]. Its special and limiting cases were discussed in Refs. [6], [53], [61], |100] . 
|102] . |144j for the simple harmonic oscillator and in Refs. [3], [12], [60], |108] . |123j for the particle 
in a constant external field; see also references therein. The damped oscillations have been studied 
to a great extent in classical mechanics [1], [5] and [79]. Their quantum analogs are introduced 
and analyzed from different viewpoints by many authors; see, for example, [13], [H], [23], [21], [25], 
[27], [32], [33], [35], [37], [10], [^, [S3], [HS], [M], [110], [130], [HO], [113], [HO], and references 
therein. The quantum parametric oscillator with variable frequency is also largely studied in view 
of its physical importance; see, for example, [22], [30], [02], [03], [SI], [05], [OZ], [117], [IIS], [121], 
|122] . jl30] . and |133] : a detailed bibliography is given in 



QUANTUM INTEGRALS OF MOTION 



3 



In the present paper we revisit a familiar topic of the quantum integrals of motion for the time- 
dependent Schrodinger equation 

^^ = H{t)iJ (1.10) 
with variable quadratic Hamiltonians of the form 

H = a(t)p^ + h it) x^ + d (t) {px + xp) , (1.11) 

where p = —id/dx, H = 1 and a{t) , b {t) , c{t) = 2d {t) are some real- valued functions of time only 
(see, for example, [SB], [HI], [HI], [HZ], [SS], |149] . |15U] and references therein). A related energy 
operator E is defined in a traditional way as a quadratic in p and x operator that has constant 
expectation values [39] : 

d d f°° 

IE) = - dx = 0. (1.12) 

dt dt 

It is well-known that such quadratic invariants are not unique. Although an elegant general solu- 
tion is known, say, for the parametric oscillator, it involves an integration of nonlinear Ermakov's 
equation |91]. Here the simplest energy operators are constructed for several integrable models of 
the damped and modified quantum oscillators. Then an extension of the familiar Lewis-Riesenfeld 
quadratic invariant is given to the most general case of the variable non-self-adjoint quadratic 
Hamiltonian (see also [SI], |149] . |15U] . we do not use canonical transformations and deal only with 
real- valued solutions of the corresponding generalized Ermakov system), which seems to be missing 
in the available literature and may be considered as the main result of this paper. (An attempt to 
collect relevant references is madqj.) Group-theoretical aspects will be discussed elsewhere, we only 
provide the factorization of the general quadratic invariant (see also |138j ). 

In general the average (E) is not positive. A complete dynamics of the expectation values of some 
energy-related positive operators is found instead for each model, which is a somewhat interesting 
result on its own. In addition to other works [7], [39], [36], [57], [91], [95], [98], [HI], [150] these 
advances allow us to discuss uniqueness of the corresponding Cauchy initial value problem for 
the special models and for the general quadratic Hamiltonian under consideration as a modest 
contribution to this well-developed area of quantum mechanics and partial differential equations. 
Further relations of the quadratic invariants with the solution of the initial value problem are 
discussed in the forthcoming paper |138] . 

The paper is organized as follows. In Section 2 we review several exactly solvable models of the 
damped and generalized oscillators in quantum mechanics. Some of these "exotic" oscillators with 
variable quadratic Hamiltonians appear to be missing, and/or are just recently introduced, in the 
available literature. The corresponding Green functions are found in terms of elementary functions. 
The dynamical invariants and quadratic energy-related operators are discussed in Sections 3 and 
4. The last section is concerned with an application to the Cauchy initial value problems. The 
classical equations of motion for the expectation values of the position operator for the quantum 
oscillators under consideration are derived in Appendix A. The Heisenberg uncertainty relation 
and linear dynamic invariants are revisited, respectively, in Appendices B and C. Solutions of a 
required differential equation are given in Appendix D to make our presentation is as self-contained 
as possible. 

complete bibliography on classical and quantum generalized harmonic oscillators, their invariants, group- 
theoretical methods and applications is very extensive. Only case of the damped oscillators in [33 includes about 
600 references! 
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2. Some Integrable Quadratig Hamiltonians 



Quantum systems with the Hamiltonians (11.111) are called the generalized harmonic oscillators 
[7], [36], [57], [M], |149] . |150j . In this paper we concentrate, among others, on the following variable 
Hamiltonians: the Caldirola-Kanai Hamiltonian of the quantum damped oscillator [13], [33], [69] , 
[146] and some of its natural modifications, a modified oscillator introduced by Meiler, Cordero- 
Soto and Suslov [lUl] , [22] , the quantum damped oscillator of Chruscihski and Jurkowski [25] in the 
coordinate and momentum representations and a quantum-modified parametric oscillator which is 
believed to be new. The Green functions are derived in a united way. 



2.1. The Caldirola-Kanai Hamiltonian. A model of the quantum damped oscillator with a 
variable Hamiltonian of the form 

H='^ (e-2A* p2 ^ ^2Xt ^2^) ^2.1) 

is called the Caldirola-Kanai model [1], [13], [33], [69], |146j . Nowadays it is a standard way of 
adding friction to the quantum harmonic oscillator. The Green function is given by 



where 



G (x, y,t) = J g<.W.H/3Wx,+,(i),^) ^ ^ ^ I 2 _ y2 > (2.2) 

\' 2TiiojQ sin tut v 
uj COS ujt — \ sin ujt 



a{t) = ^ir:^r^_^_^e^^*, (2.3) 
Zujq sm oot 

P{t) = ^e'\ (2.4) 

ooo sm oot 

UJ COS ujt + \ sin ut 

lif) = ^ ■■ 7 • (2-5) 

Zooq sm oot 

This popular model had been studied in detail by many authors from different viewpoints; see, for 
example, [2], [H], [B], [IB], [I?], [20], [21], [SZ], [ZO], [H], [72], [73], [77], [S2], [nU], [112], [115], 
|127] . |139] . |140] . |145] . |151] and references therein, a detailed bibliography can be found in [33] . 

[US]. 



2.2. A Modified Caldirola-Kanai Hamiltonian. In this paper, we would like to consider an- 
other version of the quantum damped oscillator with variable Hamiltonian of the form 



^0 f„-2Xt J2 

The Green functions in (12. 2p has 



(e-2A* p2 ^ ^2Xt ^2^) _x^px + xp) . (2.6) 



= c^cosct + Asinct^,, 
^ ^ 2cc;osincc;i ^ ^ 

(3it) = ^e^*, (2.8) 

ooo sm 

a; cos — A sin 

7 t = . 2.9 

This can be derived directly from equations (ll.2p - (ll.8p following Refs. [26] and 
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The Ehrenfest theorem for both Caldirola-Kanai models has the same form 

^ (x) + 2A^ (x) + ^x) = 0, (2.10) 

which coincides with the classical equation of motion for a damped oscillator [5], [79]. Details are 
provided in Appendix A. 

2.3. The United Model. The following non-self-adjoint Hamiltonian: 

if = ^ (e-2At p2 ^ ^2Xt ^2^ _ (2.11) 

coincides with the original Caldirola-Kanai model when /i = and the Hamiltonian is self-adjoint. 
Another special case A = corresponds to the quantum damped oscillator discussed in [2?j as an 
example of a simple quantum system with the non-self-adjoint Hamiltonian. (This is an alternative 
way to introduce dissipation of energy to the quantum harmonic oscillator.) Combining both cases 
we refer to fl2.1ip as the united Hamiltonian. 

The Green function is given by 



where 



G{x,y,t) = J ^^{ait).^+m^y^,it)y^) ^ (2.12) 

" 27iiuo sin cot 



= c^cosa;t + (/.-A)sinc.t 

2uosmut ' ^ ^ 

Ht) = ^e^*, (2.14) 

ooo sm out 

7 t = ^ ^ 2.15 



with (jj 



^ul - (A - lif > 0. 



In this case the Ehrenfest theorem takes the form: 

^ (x) + 2 (A + /.) ^ (x) + {ul + 4A/i) (x) = 0. (2.16) 
It is derived in Appendix A and the Heisenberg uncertainty relation is discussed in Appendix B. 

2.4. A Modified Oscillator. The one-dimensional Hamiltonian of a modified oscillator introduced 
by Meiler, Cordero-Soto and Suslov |101j . [28] has the form 

H = (cost J9 + sin t x)^ (2.17) 
= cos^ t + sin^ t + sin t cos t (px + xp) 

= - {p^ + x^) H — cos 2t {p^ — x^) H — sin 2t [px + xp) . 
2 2 2 
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(A physical interpretation of this Hamiltonian from the viewpoint of quantum dynamical invariants 
will be discussed in Section 4.) The Green function is given in terms of trigonometric and hyperbolic 
functions as follows 

G{x,y,t) = ^ = (2.18) 

y 2711 (cos t sinh t + sin t cosh t) 

(x^ — y"^) sin t sinh t + 2xy — (x^ + y"^) cos t cosh t 

X exp ' 



2i (cos t sinh t + sin t cosh t) 
More details can be found in |1U1] . The corresponding Ehrenfest theorem, namely, 

d 



(a;) + 2tant- (x) - 2 (x) = 0, (2.19) 



is derived in Appendix A. 



2.5. The Modified Damped Oscillator. The time-dependent Schrodinger equation 

=H{t)^P (2.20) 
with the variable quadratic Hamiltonian of the form 

H = + ^ cosh^ (At) x^ p=-^ (2.21) 

2mcosh2(At) 2 ^ ^ ' ^ 2 9x ^ ' 



has been recently considered by Chruscihski and Jurkowski [25] as a model of the quantum damped 
oscillator; see also |106j . 

In this case the characteristic equation (11. 6p takes the form 

/i" + 2A tanh (At) /i' + wg/i = 0. (2.22) 
The particular solution is given by 

h sin (wt) 



/^(t) = — ^":7::. , c^ = v^o-a'>o (2.23) 

mcj cosh (At) ^ 



and the corresponding propagator can be presented as follows 



where 



Gix,y,t) = J ^^^""^ (^^) ^i{ait).^+m^y+,(t)y^) ^ (2.24) 
y 27ri^sin [ut) 

Tfi COSll ( \t] 

a (t) = . ^ ^ (w cos (wt) cosh (At) - A sin (cot) sinh (At)) , (2.25) 
zn sin (wt j 

ma; cosh (At) 
2Aism (wt) 

, , mu cos (wt) 

7t = _ 2.27 

2^sin(cjt) ^ ^ 
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(We somewhat simphfy the original propagator found in [25]; see also |76|.) This Green function can 
be independently derived from our equations f ll.3p -f lT3]) with the help of the following elementary 
antiderivative: 

f A cos {ut + 5) sinh (At) + u sin (cut + 5) cosh (At) 
\a; cos {uot + 5) cosh (At) — A sin (cut + 5) sinh (At) 

ujuoq cosh^ (At) 
{oj cos (cut + 5) cosh (At) — A sin (cut + 5) sinh (At))^ 
Further details are left to the reader. 



(2.28) 



Special cases are as follows: when A = 0, one recovers the standard propagator for the linear 
harmonic oscillator [50], and = gives a pure damping case [76] : 



Gi^.y.t) = J , exp ( Z^^'u'/^l ] . (2.29) 

^ ^ ^' 27ri^tanh(a;t) ^ \^ 2^tanh (wt) ^ ^ ' 

In the limit A — > formula fl2.29p reproduces the propagator for a free particle [50] . 

The Ehrenfest theorem for the quantum damped oscillator of Chruscihski and Jurkowski coincides 
with our characteristic equation fl2.22p : see Appendix A for more details. 

It is worth adding that in the momentum representation, when p ^ x, a rescaled Hamiltonian 
(I2.2ip {H = muo = 1) takes the form 

(2.30) 

The corresponding characteristic equation 

/i" - 2A tanh (At) /i' + wg/i = (2.31) 
has a required elementary solution 

fj, = — (A cos (cut) sinh (At) + u sin (cut) cosh (At)) (2.32) 

(jJq 

with /i' (0) = 2a (0) = Uq and 

fi ^ —e^^ {X cos (ut) +u sin {ut)) (2.33) 
2uJo 

as t — )■ oo. The Green function is given by formula (11.20 with the following coefficients: 

^/^x ^ cjQ cos jut) 

2cosh(At) (Acos(a;t)sinh(At) +a;sin(u;t)cosh(At))' ^' ' 

pfA = ^ (2.35) 

A cos (ut) sinh (At) + u sin (cut) cosh (At) ' 

Wo (w cos (wt) cosh (At) — A sin (cut) sinh (At)) (2 36) 

2u (A cos (cut) sinh (At) + u sin (cut) cosh (At)) 

The details are left to the reader. 
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2.6. A Modified Parametric Oscillator. In a similar fashion we consider the following Hamil- 
tonian: 

H =- (tanh^ (At + 5) + coth^ (At + 6) x^) (2.37) 

+ sinh(2At + 2^) (^^ + ^^^ 
which seems to be missing in the available literature. The corresponding characteristic equation: 

4A / 2A^ A 

has an elementary solution of the form: 

. , , tanh (At + 6) 

/i = sin (ut) ^— — ^. (2.39) 

cotho 



In the limit t — )■ oo, yU — )■ sin (cut) tanh 6. 

The Green function can be found as follows 



^ y,t)-^ 27rz sin (ut) tanh (At + 5) ^ ' ^^"^"^ 

where 

a (t) = i cot (wt) coth^ (At + 6) , (2.41) 

/3(t) = -^^^coth(At + 5), (2.42) 
sm [ut) 

7 (t) = 1 cot (cot) coth^ (5. (2.43) 

The Ehrenfest theorem coincides with the characteristic equation fl2.38p . One should interchange 
a b and d — —d in the momentum representation [2B]- The corresponding solutions can be found 
with the help of the substitution 6 ^ 6 + in /2. The trigonometric cases, when A — )■ iA, 5 — ?■ z5 and 
u — )■ — are left to the reader. 



2.7. Parametric Oscillators. In conclusion a somewhat related quantum parametric oscillator: 

1/^2 f 2 2A2 

when 



H = Up'+(uj' + ) ) , (2.44) 



and 

A cos (ut) sinh (At) + u sin (wt) cosh (At) 
+ A^) cosh (At) 

has the Green function (II. 2p with the following coefficients: 



/\ L.un i^u/t; niiiii -p uj niii i^u/t; uunii 

= ^ o . ^2^ TTTT^ ' (2-46) 



(w^ + A^ cosh ^ (At)) cos (cut) — Xu tanh (At) sin (cut) 
~ 2 (w sin (ut) + A tanh (At) cos (ut)) ' ^' ^ 
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Pit) = ■ , , ^ w^ X (2.48) 

uj sm [(jjt) + A tann (At) cos [ut) 

fcj^ + A^) cos (wt) — A tanh (Ai) sin (wt)) 
^ (t) = ^ LI 1^ - 1^ - 1^ — (2 49) 

2a;(u;sin(a;t) + Atanh(At)cos(c<;t)) ' ^' ^ 

The Green function for the parametric oscillator in general: 

H =^{p^ + u\t)x^) (2.50) 

can be found, for example, in Ref. [81]. (The time-dependent quantum oscillator was thoroughly 
examined by Husimi [22] , [03] and later many authors had treated different aspects of the problem; 
see [39j, [HZ], [IS], ^T7\, [IIS], [M], [122], [ISO] and [133j; a detailed bibliography is given in 
Ref. [H].) 

3. Expectation Values of Quadratic Operators 

We start from a convenient differentiation formula. 
Lemma 1. Let 

H = a{t)p'^ + b (t) x'^ + d (t) {px + xp) , (3.1) 

O = A{t)p^ + B {t) x^ + C (t) {px + xp) (3.2) 

and 

(O) = {'ilj,Oij) = 1^*0^: dx, i-^ = Hi) (3.3) 

J — oo ^ 

(we use the star for complex conjugate). Then 



1(0) ^ (^+4(aC-.^))(.^> 

+ (f + 4(.i.-.C))(.^> 



(3.4) 



+ + 2 {aB - hA)^ {px + xp) . 
Proof. The time derivative of the expectation value can be written as [80], |102j . |129j : 

where [0.,H] = OH — HO (we freely interchange differentiation and integration throughout the 
paper, it can be justified for certain classes of solutions [22], |113] . |119] . |147] ). One should make 
use of the standard commutator properties, including familiar identities 

=2i{px + xp), [x,p'^] = ^ip, [x'^,p\ =2ix, (3.6) 
[px + xp,p'^~\ = Aip'^, [x'^,px + xp] = 4ix^, 

in order to complete the proof. □ 
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Quantum systems with the self-adjoint Hamiltonians (13. ip are called the generalized harmonic 
oscillators [7], [2S], [5Z], [HI], |149] . |15U] . At the same time one has to deal with non- self-adjoint 
Hamiltonians in the theory of dissipative quantum systems (see, for example, |27], |33], [77], |143] . 
[146] and references therein) or when using separation of variables in an accelerating frame of 
reference for a charged particle moving in an uniform time-dependent magnetic field [26]. An 
extension to the case of non-self-adjoint Hamiltonians is as follows. 



Lemma 2. // 



then 



H = a{t)p'^ + b (t) + c (t) px + d (t) xp, (3.7) 
O = A{t)p^ + B (t) x^ + C (t) px + D (t) xp, (3.8) 

d f dA 

di 



(O) = (^^ + 2a{C + D)-{3c + d)A^{p^) (3.9) 
- 26 (C + D) + (c + 3d) B^ (x^) 



Proof. One should use 



dB 

m 

+ ['^ + 2{aB -hA)-{c-d)C] {px) 
+ ['^ + 2{aB -hA) -{c-d)D] {xp) 



iO) = (^) + -{OH-H^O), (3.10) 



dt' ' \ dt 



where if^ is the Hermitian adjoint of the Hamiltonian operator H. Our formula is a simple extension 
of the well-known expression [80], |102j . |129] to the case of a non-self-adjoint Hamiltonian [27] . 
Standard commutator evaluations complete the proof. □ 



Polynomial operators of the higher orders in x and p can be differentiated in a similar fashion. 
An analog of the product rule is given in |138] . The details are left to the reader. 



4. Energy Operators and Quadratic Invariants 
In the case of the time-independent Hamiltonian, one gets 

|W=0 (4.1) 

by (13.51) . The law of conservation of energy states that 

E = {H) = constant. (4.2) 

In general one has to construct quantum integrals of motion, or dynamical invariants, that are 
different from the variable Hamiltonian (see, for example, [91], |149] . |150] : linear case is dealt with 
in [37], [39], [95], [98] and Appendix C). 



d (t) A) 


= 0, 


(4.5) 


b it) C) 


= 0, 


(4.6) 


b it) A) 


= 0. 


(4.7) 
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4.1. Energy Operators. A familiar definition is in order (see, for example, [39], [95]). 

Definition 1. We call the quadratic operator ( IJ.i^p an energy operator E , or a quadratic (dynamical) 
invariant, if 

|W=0 (4.3) 

for the corresponding variable Hamiltonian liS. 

By Lemma 1 the coefficients of an energy operator, 

E = A{t)p^ + B (t) x^ + C (t) {px + xp) , (4.4) 
must satisfy the system of ordinary differential equations: 

dA . / \ ^ 
- + 4(„WC 

dB , , 
-+4(<iWf 

In general a unique solution of this system with respect to arbitrary initial conditions Aq = A (0) , 
Bq = B [0) , Cq = C (0) [59] determines a three-parameter family of the quadratic invariants (14. 4p . 
Special cases, when solutions can be found explicitly, are of the most practical importance. 

In this section we find the simplest energy operators for all quadratic models under consideration 
as follows: 

E=^ (e-2A* p2 ^ ^2Xt ^2) + ^ (p^; + ^ (4.8) 
E = ^ (e-2At p2 + ^2Xt ^2) _ ^ + ^p) ^ (4.9) 

£' = icos2t (p^ - x^) + i sin 2t {px + px) , (4.10) 

E = tanh^ (At + S) p^ + coth^ (At + 5) x^ (4.11) 

for the Caldirola-Kanai Hamiltonian (12. ip [139] , the modified Caldirola-Kanai Hamiltonian (12. 6p . 
the modified oscillator of Meiler, Cordero-Soto and Suslov fl2.17p and for the modified parametric 
oscillator (I2.37p . respectively. Their coefficients solve the corresponding systems (14. 5p - (14. 70 for 
special initial data. 

An energy operator for the united model (12. lip is given by 

E = ^^e'^* p" + e^^* x^) + ^ (A - /x) e^* {px + xp) . (4.12) 

One should use Lemma 2; verification is left to the reader. Finally an energy operator for the 
quantum damped oscillator of Chruscihski and Jurkowski with a rescaled Hamiltonian (I5.39P is 
given by expression (I5.40p . A general case of the variable quadratic Hamiltonian is discussed in 
Theorem 1. 
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4.2. The Lewis— Riesenfeld Invariant. Classical Hamiltonian of the generalized harmonic oscil- 
lator can be transformed into the Hamiltonian of a parametric oscillator [7], [SZ], |114] . |150] . All 
quadratic invariants of the quantum parametric oscillator f l2.50p can be found as follows [HH], [89] , 
[90] , [91] . The corresponding system, 

A' + 2C = 0, (4.13) 
B' -2uj^t)C = 0, (4.14) 
C + B-uj^{t)A = 0, (4.15) 

is integrated by the substitution A = k^. Then C = —kk,', B = kk" + (k')^ + (t) and equation 
fl4.14p becomes 

{kk" + (k')^ + (t) ' + 2u;^ (t) = 0, 

h{k" + u? (t) k) ' + 3/?' (k" + u'^ it) k) = 
or with an integrating factor: 

^ [k" [k" + u"" {t) k)) = (4.16) 



(see [91] and [86]). Thus 



Co 



K" + Uj\t)K=^ (co = 0,l) (4.17) 



K 



and a general solution of the system fl4.13l) - fl4.15l) is given by 



^2 , Co 



A = f^\ 5 = (k') C = -nn' (4.18) 



in terms of solutions of the nonlinear equation (I4.17p . which is called Ermakov's equation, when 
Co = 1 [S] (see also, [85], [90], [120] and [130] ). Thus the quadratic integrals of motion can be 
presented in the form [91[ : 

E=(kp- Kxf + (4.19) 

for any given solution of the Ermakov equation fl4.17p . This quantum invariant is an analog of the 
Ermakov-Lewis integral of motion for the classical parametric oscillator [H], [88], [89], [90], [141] . 

In general if two linearly independent solutions of the classical parametric oscillator equation are 
available: 

u" + (t) u = 0, v" + {t) w = 0, (4.20) 
then solutions of the nonlinear Ermakov equation: 

1 



k" + u^{t)K = ^ (4.21) 
are given by 

K = {Au^ + 2Buv + Cv^) (4.22) 

(so-called Pinney's solution [120] . [I3] . [85] . [90] . [114] ). where the constants A, B and C are related 
according to AC—B^ = 1 /W^ with W being the constant Wronskian of the two linearly independent 
solutions. 

For example, in the case of the simple harmonic oscillator with u (t) = 1, there are two elementary 
solutions: 

K = l (co = l), ft = cost (co = 0) (4.23) 
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and the energy operators are given by 

H = ^{p' + x'), (4.24) 

E = {cost p + sint xf . (4.25) 

It provides a somewhat better understanding of the nature of the Hamiltonian discussed by Meiler, 
Cordero-Soto and Suslov |101] — this operator plays a role of the simplest time-dependent quadratic 
integral of motion for the linear harmonic oscillator. 

In a similar fashion the dynamical invariants of the parametric oscillator f l2.44p are given by the 
expression f l4.19p with Cq 7^ 0. In the Pinney solution (14.221) one can choose 

u cos (cut) cosh (At) — A sin (ut) sinh (At) , , 

Lo sin (ut) cosh (At) + A cos (ut) sinh (At) , , 

V = ^ — , , — ^ — (4.27) 

cosh (At) ^ ^ 

as two linearly independent solutions of the classical equation of motion (I2.45P with W {u, v) = 
u + A^) . If A = C and 5 = 0, then 

K= {u^ + A' tanh^ (At)) (4.28) 
is a particular solution of the corresponding Ermakov equation: 

The simplest positive energy integral for our parametric oscillator (I2.44p is given by 

E = (co' + \'t^nh'{\t)) p2 ^ ^3 sinh (At) ^ ^4_3q^ 

cosh (At) 



A^ sinh^ (At) + 00^ (A' + 00^? cosh^ (At) „ 

X . 



cosh^(At) {co^ + XHanh^ (Xt)) 
Another possibility is to take a general solution of ( 12.45^ with cq = 0. 

4.3. An Extension to General Quadratic Hamiltonians. We consider the following general- 
ization of the Lewis-Riesenfeld invariant (I4.19P (see also [51], |150] ). 



Theorem 1. The dynamical invariants for the general quadratic Hamiltonian 7| ) are given by 

E= — iKp-——x) 4.31 
jj^i \ la dt J jj2K 

where Cq is a constant, 



Hi = exp ^— J (3c + d) ds^ , /ig = exp (c + 3d) ds^ , (4.32) 



and K satisfies the auxiliary nonlinear equation: 



k^(k^] +4abk\ = ^, (4.33) 
dt \ dt J 
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where ^ 

k = ^exp (^2 (c + d) ds^ . (4.34) 

(For the self-adjoint Hamiltonians c = d.) 

The 1/2, 6 = (t) /2 and c = c? = 0, corresponds to the original invariant fl4.19p . 

Proof. By Lemma 2 in order to find quadratic invariants of the form 

E = Ap^ + Bx^ + Cpx + Dxp (4.35) 
we have to solve the following system of ordinary differential equations: 



d A 

- + 2a(C + Z,)- 


{3c + d)A = 


0, 


(4.36) 




{c + 3d)B = 


0, 


(4.37) 


dC 

— + 2(aB- bA) - 
dt 


-{c-d)C = 


0, 


(4.38) 


^ + 2(aB- bA) - 
dt 


-{c-d)D = 


0, 


(4.39) 


say, for arbitrary analytic coefficients a{t) , b (t) , 


c (t) and d (t) . 


The substitution C 





D = Ci — Di allows one to transform the last two equations: 

dC 

—^ + 2{aB- bA) -{c-d)Ci = 0, (4.40) 
dt 

= {c- d) Di, Di = constant exp {c - d) ds^ . (4.41) 

Then 

Cpx + Dxp = Ci {px + xp) + Di {px — xp) 

and, in view of the canonical commutation relation, the coefficient Di can be eliminated from the 
consideration as belonging to the linear invariants (see appendix C). 

Introducing integrating factors into f l4.36p . f l4.37p and f l4.40p . we get 

d u' 

— {^i^A) + 4a/iiCi = 0, ^ = -3c - d, (4.42) 
dt 

^ {fi^B) - Abfi.Ci = 0, ^ = c + 3rf, (4.43) 
dt jj,2 

4- (yUgCi) + 2/^3 {aB - M) = 0, i^ = -c + d (4.44) 
dt /i3 

with /i| = yUi/i2- After the substitution 

A = ix,A, B = ^2B, C = ix.,Ci, (4.45) 

the system takes the form 

^ + 4aMc = 0, (4.46) 
dt y fj,2 
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''^ +2(aMB-bMA]=0. (4.48) 



Introducing a "proper time" : 



we finally obtain: 



dt \ y (j'2 V /^i 



T= I 2aJ^^ds, (4.49) 



V + 2C = 0, (4.50) 
dr 

dR ~ 

^-2uj''{t)C = (4.51) 
dr 

^ + 5_^2(^)2=0, uj\t)=^-^, (4.52) 

dr ajJi 

which is identical to the original Lewis-Riesenfeld system (I4.13p - fl4.15l) (positivity of cj"^ is not 
required). The solution is given by 

where k satisfies the Ermakov equation: 

'^..Hr)..% .^M.^, (4.54) 
with respect to the new time fl4.49p . In view of 



d ^d 
dr dV 



k = ^exp (2 (c + d) ds^ , (4.55) 



the Ermakov equation fl4.54p is transformed into our auxiliary equation fl4.33p . A back substitution 
results in the dynamical invariant f l4.3ip when the square is completed. □ 

Lemma 3. The dynamical invariant ( [^ . 31\ ) can be represented in more symmetric form 

dfi i\ \ \ ^0 2 



X exp yj {c — d) ds 
where Co is a constant and fi is a solution of the following auxiliary equation: 

Ij" --fi' + (Aab+ (--€- d] (c + d)-c'- d') 12 = Co^^. (4.57) 
a \ \a J J jji-^ 
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Proof. Use the substitution 

K = fiexp^—J (c + d) ds^ (4.58) 
in f l4.3ip and fl4.33p . A somewhat different proof is given in |138] . □ 

The corresponding classical invariant is discussed, for example, in Refs. |141] and |150] . (Compare 
also our expression ( I4.56P with the one given in the last paper for the self-adjoint case; we give a 
detailed proof for the non-self-adjoint Hamiltonians and emphasize connection with the Ermakov 
equation.) 

It is worth noting, in conclusion, that, if /i^ and /ig are two linearly independent solutions of the 
linear equation: 

/i" - -fi' + (iab +(--c-d^{c + d)-c'-d'^fi = 0, (4.59) 



the general solution of the nonlinear auxiliary equation (I4.57P is given by 

/X = {Af4 + 25/ii/i2 + Cfxl) ^''^ , (4.60) 
where the constants A, B and C are related according to 

with W (/ii,/i2) = 1^\IJ^2 ~ A*'i/^2 = constant (2a) being the Wronskian of the two linearly indepen- 
dent solutions. This is a simple extension of Finney's solution fl4.22p : our equations fl4.57p and 
f l4.59p form the generalized Ermakov system |43] , |114j . Further generalization of the superposition 
formula fl4.60p - fl4.6ip is discussed in Ref. |138j . (If Co ^ 0, the substitution jj — )■ C^^^fi reduces 
equation f l4.57p to a similar form with Cq = 1.) Special case of the time- dependent damped harmonic 
oscillator, when a = e~-^'^^'^ /2, h = u"^ {t) €^^^^2, F (t) = f (s) ds and c = c? = 0, is discussed in 
[82], [83]. 

4.4. An Example. The simplest energy operators have been already discussed in section 4.1 for 
all models of quantum oscillators under consideration. In order to demonstrate how the general 
approach works we discuss the united Hamiltonian (12. lip , when a = (u;o/2) e~^^*, b = (a;o/2)e^'^* 
and c = 0, d = — /i. A direct calculation shows that the function 




. = e-- (4.62) 



satisfies the following equation 



k" + 2Xk' + uIk=('^Y^, = ul-{\- fif >0, (4.63) 
V 2 / K'^ 

which corresponds to the nonlinear auxiliary equation (14.570 with Co = u'^/A. The quadratic invari- 
ant fl4.56p simplifies to the previously found expression fl4.12p . Solution f l4.60p can be used for the 
most general case. Details are left to the reader. 
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4.5. Factorization of the Dynamical Invariant. Following Ref. |27] the energy operator fl4.56p 
can be presented in the standard harmonic oscillator form: 

E ='^[a{t)aHt) + a^t)a{t)) , (4.64) 



where 



u {t) = uoexp (c - d) ds^ , uq = 2a/Co > 0, (4.65) 
\ 2fi 2a^fijjQ ) ^/(Uo^x' 



2yU 2a^/^Jo J y/uJo dx 

and /i is a solution of the nonlinear auxiliary equation ( I4.57P . Here the time-dependent annihilation 
a (t) and creation (t) operators satisfy the usual commutation relation: 

a{t)a^ (t) -a^t)a{t) = 1. (4.68) 

The oscillator-type spectrum and the corresponding time- dependent eigenfunctions of the dynamical 
invariant E can be obtain now in a standard way by using the Heisenberg-Weyl algebra of the rasing 
and lowering operators (a "second quantization" [91], the Fock states). Explicit solution of the 
Cauchy initial value problem in terms of the quadratic invariant eigenfunction expansion is found 
in Ref. |138] . In addition the n-dimensional oscillator wave functions form a basis of the irreducible 
unitary representation of the Lie algebra of the noncompact group SU (1,1) corresponding to the 
discrete positive series P;^ (see |101] . [Ill] and |132] ). Our operators f l4.66p -( l4.67p allow one to 

extend these group-theoretical properties to the general dynamical invariant fl4.64p . We shall further 
elaborate on these connections elsewhere. 

5. Application to the Cauchy Initial Value Problems 

Explicit solution of the initial value problem in terms of eigenfunctions of the general quadratic 
invariant is given in Ref. [138] . Here we formulate the following uniqueness result. 

Lemma 4. Suppose that the expectation value 

{Ho) = {iJ,Hoij)>0 (5.1) 

for a positive quadratic operator 

Ho = f (t) (a {t)p + f3 (t) xY + g (t) x' (/ (t) >0, g (t) > 0) (5.2) 
(a (t) and /3 {t) are real-valued functions) vanishes for all t G [0, T) : 

{Ho) = (Ho) (t) = (Ho) (0) = 0, (5.3) 
when ip {x, 0) = almost everywhere. Then the corresponding Cauchy initial value problem 

'% = Hip, iP{x,0) = <f{x) (5.4) 

may have only one solution, when xip {x,t) G (M) (if Hq = g{t) /, where I = id is the identity 
operator, ip & L"^ (R) ). 
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Here it is not assumed that Hq is the quantum integral of motion when j| (Hq) = 0. 



Proof. If there are two solutions: 



dt dt 

with the same initial condition ipi {x, 0) = ip2 i^y 0) = ? then by the superposition principle the 
function ip = ip-^—ilj2 also a solution with respect to the zero initial data ip {x, 0) = ip {x)—ip {x) = 
0. By the hypothesis of the lemma 

{tP, Hoij) = f (t) {{ap + /3x) ij, {ap + (3x) t/j) + g (t) {x^j, xt/j) = 

for all t G [0,T). Therefore xip {x,t) = x {ipi {x,t) — ip2 {x,t)) = and ipi {x,t) = i^yt) almost 
everywhere for alH > by the axiom of the inner product in (M) . □ 

In order to apply this lemma to the variable Hamiltonians one has to identify the corresponding 
positive operators Hq and establish their required uniqueness dynamics properties with respect to 
the zero initial data. In addition to the simplest available dynamical invariant (IB. II) . it is worth 
exploring other (quadratic) possibilities. The authors believe that it is interesting and may be 
important on its own. For example, our approach gives an opportunity to determine a complete 
time-evolution of the standard deviations (]B.9|) - (]B.10p for each of the generalized harmonic oscil- 



lators under consideration. The details will be discussed elsewhere. 

5.1. The Caldirola-Kanai Hamiltonian. The required operators are given by 

H = Ho='^ p^ + e^^* x^) , (5.5) 

L = ^ = Xujo i-e-^'' p' + e''' x') , (5.6) 
E = ^ p^ + e'^' x^) + ^ipx + xp) , j^{E) = 0. (5.7) 



By ([33]) 



Applying formula ( 13. 4 p one gets 

d 



and 



^^{L) = 2X'coo{e'''' {p')+e''' {x')) (5.9) 
+2A(X'o {px + xp) 

^^{L)+4co'{H)=Au;l{E), (5.10) 



with the help of ([531) and (ISTj) . 

In view of (15.81) and (I5.10p the dynamics of the Hamiltonian expectation value (H) is governed 
by the following second-order differential equation 



^^{H)+4uj'{H)=4ul{E), (5.11) 
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with the unique solution given by 

(^jj^ ^ \ h^^ o\ /o (2a;t) + _ ^ — ^^ sin {2u;t) + -f {E), . (5.12) 

The hypotheses of Lemma 4 are satisfied. Our solution allows to determine a complete time- 
evolution of the expectation values of the operators p^, and px + xp. Further details are left to 
the reader. 



5.2. The Modified Caldirola-Kanai Hamiltonian. The required operators are 

= ^ p2 + e^"* x') -X{px + xp) , (5.13) 



^ = Ac.oH--p^ + e-.^) = ^, (5.14) 



E = ^ (e-2At p2 ^ g2At ^2) _ ^ (pa; + a;p) . (5.15) 



We consider the expectation value (Hq) of the positive operator 



H, = ^ {e-''' p' + e''' x') . (5.16) 



In this case H = 2E — Hq, and 



|{L> = 4u^H,}~4u,l{E)„, (5.18) 



which results in the differential equation (15. lip with the explicit solution 



tu^{Ho),-tul{E)^ 1 /dHo\ ...... u;, 





of the initial value problem. The hypotheses of the lemma are satisfied. 

5.3. The United Model. The related operators can be conveniently extended as follows 

Ho = ^e'^* p2 ^ g2At ^2^ ^ (5 20) 

M = e''* {px + xp) (5.22) 

and 

E = Ho{t) + ^{X-fi)M{t) (5.23) 

Then by Lemma 2 



i^gM (e-2Ai p2 ^ g2Ai ^2) + ^ (A - ^) e^* (px + xp) . 



^ (M) = -2uo (L) , (5.24) 
at 

I (//o) = (A - ^) (L) , (5.25) 
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I (S) = (5.26) 

and 

^ (L) = 4^ {Ho) + 2a;o (M) . (5.27) 

In terms of the energy operator 
and as a result 

^{Ho)+Au^Ho)^AujI{E),, u^^ujl-{\-iif>Q (5.29) 
with the unique solution of the initial value problem given by 



4(A-A^)^;^(L)oSin(2a;i) + ^(i^)o. 
The hypotheses of Lemma 4 are satisfied. 



where 



5.4. The Modified Oscillator. The required operators are 

H — {cost p + sin t x)"^ (5.31) 

= cos^ t + sin^ t + sin i cos t {px + xp) 

11 1 
= - (p^ + x^) H — cos2t (p^ — a;^) H — sin2t (px + px) 

= Ho + E{t), 

1 1 

E^E{t) = - cos 2t (p2 -x"^) +- sin 2t (pa; + px) (5.33) 

dE 

L = — = — = - sin 2t (p^ - x^) + cos 2t (px + px) . (5.34) 
at at 

I (L) = 4 {Ho) . (5.36) 



and 



Here 



and 



The expectation value {Ho) satisfies the following differential equation 

rf2 



with the explicit solution 



{Ho) = 4 {Ho) (5.37) 



{Ho) = (//o)o cosh (2t) + ^ (L)o sinh {2t) . (5.38) 
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The hypotheses of Lemma 4 are satisfied. 

5.5. The Modified Damped Oscillator. Let H = muo = 1 in the Hamiltonian (12.211) : 
without loss of generality. The corresponding energy operator can be found as follows 

E = "° y + + (540) 

2 cosh' ( At) ^ 2uo ^ ^ 

+ — tanh (At) {px + xp) , — {E) = 0, 

^ (JjV 

in view of (I4.5p - (l4.7p (one should replace A -H- 5, C — )• —C in the momentum representation). 
Introducing the following complementary operators 

2 

cosh (At) 

L = ^ cosh2(At) x^, (5.42) 

cosh2(At) V ; , V ; 

M = px + xp, (5.43) 



we get 



Then 



^{Ho) = -2Atanh(At) (L) , (5.44) 

^ (L) = -2A tanh (At) {Hq) - 2uo (M) , (5.45) 

^(M) = 2uo{L). (5.46) 

E = ^(l- ^ ^\ ]ho + L (5.47) 

2 V 2^2 posh' (At) y 4wocosh'(At) ^ ^ 

+ ^ tanh (At) M 

and, eliminating (M) and (L) from the system, one gets: 

^ (^o) - -T^^ (Ho) + 2 (2u' + —i—-) {Ho) = SCO {E), . (5.48) 



(L). (5.49) 



dt^ ' sinh (2At) dt' \ cosh' (At) 

The required initial conditions: 

d 



dt 



{Ho)^ =0, (^coth(At)^(ifo)) =-2A 



follow from (I5.44p . The unique explicit solution is given by 

. -A ^'<^>;!-°fff° (5.50) 
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(2a;tanh (At) sin {2ujt) + A (l + tanh^ (At)) cos {2ut)) 



+2 {E), 



^0 



V 2a;g cosh^ (At) 
(see appendix D). The hypotheses of Lemma 4 are satisfied. 

5.6. The Modified Parametric Oscillator. In the case fl2.37p . the energy operator (14. lip is a 
positive operator: 

{E) = tanh^ (At + 6) {p^) + coth^ (At + 5) (x^) = (E)^ > 0. (5.51) 
The related operators are 

L = tanh^ (At + 5) p2 _ coth^ (At + 5) x^ (5.52) 
M = px + xp, (5.53) 



with 

From here 



I (L) = -2u (M) , I (M) = -2u (L) . (5.55) 

^ (L) + Au' {L) = 0, ^ (M) + 4C.2 (M) = 0, (5.56) 
which determines the time-evolution of the expectation values. 

5.7. Parametric Oscillators. In general the Lewis-Riesenfeld quadratic invariant (14.191) for the 
parametric oscillator f l2.50p is obviously a positive operator for real-valued solutions of the Ermakov 
equation fl4.17p that satisfies the conditions of our lemma. 



5.8. General Quadratic Hamiltonian. In the case of Hamiltonian (13. 7p applying formula (13. 9p 
to the operators, O = {p"^, x'^,px + xp} , one obtains [27] : 

/ ip') \ /-3c(t)-rf(t) -26 (t) \ / ip') \ 

^ = c (t) + 3rf (t) 2a (t) . (5.57) 

V {px + xp) ) \ 4a (t) -46 (t) -c (t) + d{t) ) \ {px + xp) ) 

This system has a unique solution for suitable coefficients [SH] , which allows one to apply Lemma 4, 
say, for the positive operator x^. Our Theorem 1 provides another choice of positive operators. On 
the second thought a positive integral (IB.3P determines time-evolution of the squared norm and 



guarantees uniqueness in (R) . Details are left to the reader. 
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Appendix A. The Ehrenfest Theorems 

Application of formula ( 13. 5 p to the position x and momentum p operators allows one to derive 
the Ehrenfest theorem [12], |102] . |129] for the models of oscillators under consideration. 

For the Caldirola-Kanai Hamiltonian (12. ip one gets 

^ {x) = cooe-'"' ip) , I (p) = -uoe''' {x) . (A.l) 

Elimination of the expectation value (p) from this system results in the classical equation of motion 
for a damped oscillator [5], j75] : 

^ (x) + 2a| (x) + ul (x) = 0. (A.2) 
For the modified Caldirola-Kanai Hamiltonian (12. 60 the system 

j- {x) = cooe-'^' (p) - 2A (x) , ^(p) = -cooe''' {x) + 2A (p) (A.3) 
at at 

gives the same classical equation. 

In the case of the united model (12.110 one should use the differentiation formula (I3.10p . Then 

I (x) = cooe-''' ip) - 2/i (x) , j^{p) = -cooe''' (x) (A.4) 
and the second order equations are given by 

^{x)+ 2{\ + ii)j^ (x) + [ul + 4A^) (x) = 0, (A.5) 

^ (p) + 2 - A) I {p) + ujI {p) = 0. (A.6) 

The general solutions are 

(x) = Ae~(^+^)* sin (wt + 5) , (A.7) 

(p) = Ee^^-'^)* sin (wt + 7) , (A.8) 



where oj = \J — {\ — jj) > 0. 

In a similar fashion for a modified oscillator with the Hamiltonian (I2.17P we obtain 

^(x) = 2cos^t (p) + 2sin)f:cost (x) , (A. 9) 
at 

^{p) = -2sin^t (x) - 2sintcost (p) . (A.IO) 

Then 

d'^ d 

— (x) + 2tant— (x) - 2 (x) = 0, (A.ll) 

dt^^ ' dt^ ' ^ ' ^ ' 

which coincides with the characteristic equation (II. 6p in this case [28]. 
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In the case of the damped oscillator of Chruscihski and Jurkowski one obtains 



1(3.) = (A 

dt mcosh^ (At) 

^{p) = -mul cosii^ (Xt) {x) . (A.13) 

The Ehrenfest theorems coincide with the Newtonian equations of motion [25] : 

(P d 

— (x) + 2A tanh (At) — (x) + ul (x) = 0, (A.14) 

(Jjv LLL 

^ (p) - 2A tanh (At) | {p) + ul {p) = (A.15) 
with the general solutions given by 



{p) = 5 (A COS (cut + 5) sinh (At) + w sin (wt + 5) cosh (At)) , (A. 17) 



respectively. It is worth noting that both equations (IA.2P and flA.14p give the same frequency of 
oscillations for the damped motion; see [25] for more details. 

Combining all models together for the general quadratic Hamiltonian (\3.7\\ : 

I (x) = 2a it) (p) + 2d (t) (x) , I (p) = -2b (t) (x) - 2c (t) (p) (A. 18) 

with the help of fl3.10p . The Newtonian-type equation of motion for the expectation values has the 
form 

^(x)-r(t)|(x)+4a(t)(x) = (A.19) 



with 



r (t) = - - 2c + 2rf, a(t) = ab-cd + - (--^] . (A.20) 
a 2 \ a d 



In order to explain a connection with the characteristic equation fll.6l) - fll.7l) we temporarily replace 
c — )■ Co and d ^ d^ in the original Hamiltonian (11. ip . Then it takes the standard form (13.71) . if 
Cq = c + d and d^ = c. Using the new notations in fll.6l) - fll.7l) we find 

r-ro = 4(rf-c), a - ao = ^ (^^J . (A.21) 

Therefore our characteristic equation (11. 6p coincides with the corresponding Ehrenfest theorem 
( IA.19P only in the case of self-adjoint Hamiltonians, when c = d (or cq = 2do)- The united model 
shows that these equations are different otherwise. 
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Appendix B. The Heisenberg Uncertainty Relation Revisited 

A detailed review with an extensive list of references is given in Refs. f38] and [109] (see also [99]). 
We only discuss the Heisenberg uncertainty relation for the position x and momentum p = —id/dx 
operators (in the units of H) in the case of the general quadratic Hamiltonian ( 13 .Tp . By our Lemma 2 
the simplest integral of motion is given by 

Eq = exp (c(r) — (i(r)) dr^ {px — xp) (B.l) 

with 

[x,p\ = xp — px = i. (B.2) 
This implies the following time evolution: 

^) = exp ( f {d (r) - c (r)) dr] (B.3) 



of the squared norm of the wave functions. 
With the expectation values 



and the operators 
let us consider 



Ax = x — x, Ap = p — p (B.5) 



< {{Ax + iXAp)^,{Ax + iXAp)ip) (B.6) 
= {ijj, {Ax - iXAp) {Ax + iXAp) ip) 
= {{Axf)-X{l) + X'{{Apf) 
for a real parameter A. Here we have used the operator identity 

(Ax - iXAp) {Ax + iXAp) = {Axf - A + {Apf . (B.7) 

Then one gets 



t 



[{Apf) {{Axf) > i (1)^ = ^exp I 2 / (d(r) - c(r)) dr | , (B.8) 



if (1)q = {iP,4')q = 1. For the standard deviations: 



.2 {(^Pf) 



(Spy = = {P') - {pY , (B.9) 

{5xf = ^^^^=W)-{x)\ (RIO) 

we finally obtain 

6p6x>^ (B.ll) 

in the units of H. It is worth noting that Eq. fIB.lip is derived, in fact, for any operators x and 
p with the commutator ( ]B.2p — the structure of the quadratic Hamiltonian (13 .yp has only been 
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used in the norm (IB.Sp . Time-evolution of the standard deviations flB.9p - flB.10p will be discussed 
elsewhere. 

Appendix C. Linear Integrals of Motion: The 

Dodonov-Malkin-Man'ko-Trifonov Invariants 

All invariants of the form 

P = A{t)p + B{t)x + C{t) (C.l) 

for the general quadratic Hamiltonian (13. 7p can be found as follows (see, for example, [37], [39] . 
[95] . [98] and references therein). Use of the differentiation formula fl3.10p results in the following 
system: 

^ = 2c(t)A-2a(t)B, (C.2) 
at 

^ = 2b{t)A-2d{t)B, (C.3) 

^ = {cit)-dit))C. (C.4) 

The last equation is explicitly integrated and elimination of B and A from flC.2p and f lC.3p . respec- 
tively, gives the second-order equations: 

A''-(^^ + 2c~2c?jA' + A(^ab-cd+^(^^-^^^A = 0, (C.5) 

B" - (^^ + 2c-2d^B' + 4(^ab-cd-^(^^~^^^B = 0. (C.6) 

The first is equivalent to our characteristic equation f ll.6p -f ir7rp and coincides with the Ehrenfest 
theorem flA.19p - flA.20p when c ^ d. 

Thus the linear quantum invariants are given by 

P = A{t)p+ ^'^^^^^^^I'^^^^^ a: + Co exp (c (r) - d (r)) dr^ , (C.7) 

where A (t) is a general solution of equation flC.SP depending upon two parameters and Cq is the 
third constant. Our Theorem 1 gives a similar description of the quadratic invariants in terms of 
solutions of the auxiliary equation f l4.33p . Relations between linear and quadratic invariants are 
analyzed in |138] . Group-theoretical applications are discussed in [22], [36], [39], [95], [37], [91], [98] 
and elsewhere. 

Appendix D. An Elementary Differential Equation 

The nonhomogeneous differential equation of the form 

4A , / o 2A 



- sinh(2At + 27)^' + T ^ cosh^ (A^ + 7) J ' = ' ^""''^ 
{u, A and 7 are some parameters) has the following general solution: 

y = Cm{t) + C2y2{t) + Y{t), (D.2) 
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where Ci and C2 are constants, 

Ui = uj tanh (At + 7) cos {uot) — A (l + tanh^ (At + 7)) sin (cut) , (D.3) 

y2 = oj tanh (At + 7) sin {ut) + A (l + tanh^ (At + 7)) cos (wt) (D.4) 

are the fundamental solutions of the corresponding homogeneous equation with the Wronskian given 
by 

W {yi, 1/2) =uj{u^ + 4A') tanh^ (At + 7) , (D.5) 

and 

if 2A^ \ 

^ ^ ^ 1^^ " + ^x') cosh^ (At + 7) j ^^'^^ 
is a particular solution of the nonhomogeneous equation. 

One can also verify that functions: 

Zl = cos (cut) — A coth (At + 7) sin (wt) , (D.7) 
Z2 = a; sin (cut) + A coth (At + 7) cos (cut) (D.8) 

with W (zi, Z2) = CO (cu^ + A^) are fundamental solutions of the following equation: 

z" +(cu^ -1 z = (D.9) 

V sinh2(At + 7); ^ ^ 

and then carry out the substitution y = z tanh (At + 7) . Details are left to the reader. The particular 
solution of the nonhomogeneous equation can be found by the variation of parameters and/or verified 
by the substitution. Review of other integrable second-order differential equations is given in [H] . 
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